home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2001 May / SGI IRIX Base Documentation 2001 May.iso / usr / share / catman / p_man / cat3 / complib / TINVIT.z / TINVIT
Encoding:
Text File  |  1998-10-30  |  5.4 KB  |  133 lines
  1.  
  2.  
  3.  
  4. ____TTTTIIIINNNNVVVVIIIITTTT((((3333FFFF))))                                                        ____TTTTIIIINNNNVVVVIIIITTTT((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      TINVIT, STINVIT  -  EISPACK routine.  This subroutine finds those
  10.      eigenvectors of a TRIDIAGONAL SYMMETRIC matrix corresponding to specified
  11.      eigenvalues, using inverse iteration.
  12.  
  13.  
  14. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  15.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee  ttttiiiinnnnvvvviiiitttt((((nnnnmmmm,,,,nnnn,,,,dddd,,,,eeee,,,,eeee2222,,,,mmmm,,,,wwww,,,,iiiinnnndddd,,,,zzzz,,,,iiiieeeerrrrrrrr,,,,rrrrvvvv1111,,,,rrrrvvvv2222,,,,rrrrvvvv3333,,,,rrrrvvvv4444,,,,rrrrvvvv6666))))
  16.           iiiinnnntttteeeeggggeeeerrrr           nnnnmmmm,,,, nnnn,,,, mmmm,,,, iiiieeeerrrrrrrr,,,, iiiinnnndddd((((mmmm))))
  17.           ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn dddd((((nnnn)))),,,, eeee((((nnnn)))),,,, eeee2222((((nnnn)))),,,, wwww((((mmmm)))),,,, zzzz((((nnnnmmmm,,,,mmmm))))
  18.           ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn rrrrvvvv1111((((nnnn)))),,,, rrrrvvvv2222((((nnnn)))),,,, rrrrvvvv3333((((nnnn)))),,,, rrrrvvvv4444((((nnnn)))),,,, rrrrvvvv6666((((nnnn))))
  19.  
  20.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ssssttttiiiinnnnvvvviiiitttt((((nnnnmmmm,,,,nnnn,,,,dddd,,,,eeee,,,,eeee2222,,,,mmmm,,,,wwww,,,,iiiinnnndddd,,,,zzzz,,,,iiiieeeerrrrrrrr,,,,rrrrvvvv1111,,,,rrrrvvvv2222,,,,rrrrvvvv3333,,,,rrrrvvvv4444,,,,rrrrvvvv6666))))
  21.           iiiinnnntttteeeeggggeeeerrrr           nnnnmmmm,,,, nnnn,,,, mmmm,,,, iiiieeeerrrrrrrr,,,, iiiinnnndddd((((mmmm))))
  22.           rrrreeeeaaaallll             dddd((((nnnn)))),,,, eeee((((nnnn)))),,,, eeee2222((((nnnn)))),,,, wwww((((mmmm)))),,,, zzzz((((nnnnmmmm,,,,mmmm))))
  23.           rrrreeeeaaaallll             rrrrvvvv1111((((nnnn)))),,,, rrrrvvvv2222((((nnnn)))),,,, rrrrvvvv3333((((nnnn)))),,,, rrrrvvvv4444((((nnnn)))),,,, rrrrvvvv6666((((nnnn))))
  24.  
  25.  
  26.  
  27. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  28.      On Input
  29.  
  30.      NNNNMMMM must be set to the row dimension of two-dimensional array parameters
  31.      as declared in the calling program dimension statement.
  32.  
  33.      NNNN is the order of the matrix.
  34.  
  35.      DDDD contains the diagonal elements of the input matrix.
  36.  
  37.      EEEE contains the subdiagonal elements of the input matrix in its last N-1
  38.      positions.  E(1) is arbitrary.
  39.  
  40.      EEEE2222 contains the squares of the corresponding elements of E, with zeros
  41.      corresponding to negligible elements of E.  E(I) is considered negligible
  42.      if it is not larger than the product of the relative machine precision
  43.      and the sum of the magnitudes of D(I) and D(I-1).  E2(1) must contain
  44.      0.0e0 if the eigenvalues are in ascending order, or 2.0e0 if the
  45.      eigenvalues are in descending order.  If  BISECT, TRIDIB, or  IMTQLV  has
  46.      been used to find the eigenvalues, their output E2 array is exactly what
  47.      is expected here.
  48.  
  49.      MMMM is the number of specified eigenvalues.
  50.  
  51.      WWWW CONTAINS the M eigenvalues in ascending or descending order.
  52.  
  53.      IIIINNNNDDDD contains in its first M positions the submatrix indices associated
  54.      with the corresponding eigenvalues in W -- 1 for eigenvalues belonging to
  55.      the first submatrix from the top, 2 for those belonging to the second
  56.      submatrix, etc.  On Output
  57.  
  58.      ZZZZ contains the associated set of orthonormal eigenvectors. any vector
  59.      which fails to converge is set to zero.
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. ____TTTTIIIINNNNVVVVIIIITTTT((((3333FFFF))))                                                        ____TTTTIIIINNNNVVVVIIIITTTT((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      IIIIEEEERRRRRRRR is set to Zero       for normal return, -R         if the
  75.      eigenvector corresponding to the R-th
  76.         eigenvalue fails to converge in 5 iterations.
  77.  
  78.      RRRRVVVV1111 , RV2, RV3, RV4, and RV6 are temporary storage arrays. Questions and
  79.      comments should be directed to B. S. Garbow, APPLIED MATHEMATICS
  80.      DIVISION, ARGONNE NATIONAL LABORATORY
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.